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We present the results of a search for cosmic ray point sources at energies above 1.0 x 10 19 eV in the HiRes 
stereo data set. The analysis is based on a maximum likelihood ratio test using the probability density function 
for each event rather than requiring an a priori choice of a fixed angular bin size. The search is extended to the 
combined data set of HiRes data above 1.0 x 10 19 eV and AGASA data above 4.0 x 10 19 eV. In both cases, no 
CN| ■ statistically significant clustering of events consistent with a point source is found. 
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1. Introduction 



When analyzing skymaps of cosmic ray arrival directions, it is important to correctly account for the error on 
the angular resolution of individual events. This is especially true in cases where data from detectors with 
■ different angular resolution are combined, or when the angular resolution of the detector varies strongly with 

energy or event geometry. A method that correctly accounts for the angular resolution on an event-by-event 
basis was recently introduced in 1 1 1. In this paper, we perform a likelihood ratio test of the hypothesis that 
several events in the skymap come from a common source. In other words, we test whether any given position 
on the sky harbors a source which contributes n s > 1 source events to the data set. The likelihood of this 
l/"") \ hypothesis is compared to the null hypothesis n s = and this likelihood ratio is maximized using n s as a free 

parameter. By calculating the likelihood ratio for a dense grid of points on the sky, we essentially search the 
entire sky for the most likely position of a source of n s events. The statistical significance is then estimated 
fS ^ by applying the same method to a large set of random isotropic data sets and evaluating what fraction of them 

have a likelihood ratio which is equal to or larger than the ratio observed in the real data. 

In 1 1 1, this method is applied to the combined HiRes and AGASA data set of events with energies in excess 
of 4 x 10 19 eV. No excess is found. The strongest signal, from the position of the AGASA "triplet," Q is 
consistent with random clustering even if the HiRes threshold is lowered a posteriori to include a HiRes event 
with energy 3.76 x 10 19 eV in the vicinity of the "triplet." Nevertheless, a point source at that position has 
been claimed (HID, partly based on the fact that another HiRes event with energy between 1.0 x 10 19 eV 
and 3.0 x 10 19 eV is found in the same region. Here, we extend the energy range of this analysis down to 
1.0 x 10 19 eV to understand whether the additional low energy event strengthens the case for a point source 
or is consistent with expectations for an isotropic skymap with 27 1 events, the number of events in the HiRes 
stereo data set above 1.0 x 10 19 eV. We will analyze the HiRes stereo data set alone and in combination 
with the AGASA data set, with the understanding that every "signal" found in the combined set would be 
strongly biased as the AGASA threshold of 4.0 x 10 19 eV was chosen in |5| precisely because it maximized 
the clustering signal. Also, magnetic smearing, which is likely to become an issue for energies as low as 
1.0 x 10 19 eV, is not taken into account. 



2. Analysis 



We give a short description of the likelihood method used in this analysis. Consider first a fixed source location 
in right ascension and declination: x s = (a s , S s ). Given a sample of N cosmic ray events, we suppose that n s 
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events come from the source location x s , and that the remaining N — n s events are random background events. 
If the ith event is a source event, then its true arrival direction is x s . The probability density for finding it at 
some location x is given by the function Qi(x,x s ), derived from the angular errors of the event and the angular 
displacement between x s and x. On the other hand, if the ith event is a background event, then the probability 
density for finding it at a location x is given by the function Ri(x), derived from the relative exposure of the 
detector to an isotropic background of cosmic rays. The subscript is necessary because R depends on whether 
event i is a HiRes or AGAS A event. Each of these functions is normalized to unity over positions x in the sky. 

We do not hypothesize which individual events are source or background events. We only suppose that there are 
n s events in the sample that come from some source position x s . Therefore, the partial probability distribution 
of arrival directions x for the ith event is given by: 

Tl N — ■ 71 

Pi(x,x s ) = -j^ Qi(x,x s ) H j^-Ri(x) . (1) 

It follows that the probability of finding the ith event at the location Xi (where it is actually observed) is 

Pi^Xi, X s ). 

The likelihood for the entire set of N events is then given by: 

N 

C(n s ,x s ) =Y[Pi{xi,x s ) . (2) 

The best estimate for the number of source events, under the assumption of a point source located at x s , is 
determined by finding the value of n s which maximizes C 

Because the value of the likelihood function depends on the number of events, a more useful quantity than C 
is the likelihood ratio 1Z: 
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where C(0,x s ) is the likelihood function of the null hypothesis (n s — 0). In practice, we maximize In TZ, 
which is equivalent to maximizing C. 

The method described so far is sufficient for testing a specific source position x s . To search the entire sky for 
the source position with the strongest signal, we calculate the likelihood ratio In 1Z for a dense grid of points 
on the sky covering the full range of equatorial coordinates accessible to AGASA and HiRes. The source 
position is essentially treated as a free parameter, along with the number of source events n s . Searching for the 
parameters a s , 5 S , and n s which maximize the likelihood ratio will therefore give us the best estimate for the 
position of the source and the number of events it contributes. 

The search proceeds as follows. The visible region of the sky is divided into a fine grid of points with separa- 
tions 0.1°/ cos 5 and 0.1° in a and 5, respectively. Each point is treated in turn as a source location x s . At each 
point, the specific quantities Qi(xi,x s ) and Ri (xi ) are required for every event. The source probability density 
function Qi(xi, x s ) depends on the angular resolution associated with the zth event, which in principle may 
depend on several quantities including energy, zenith angle, etc. The background probability density function 
Ri(xi) depends on the detector exposure to different parts of the sky: it is generally the same function for all 
events observed by a given detector, but may in principle be a function of e.g. energy as well. For each event, 
Qi(x"i, x s ) needs to be reevaluated at every grid point, while Ri(xi) needs only to be evaluated one time. 



Search for UHECR Point Sources in HiRes Stereo Data 



3 



Once the values of Qi and Ri are specified, the log likelihood ratio In 1Z is maximized with respect to n s , 
where n s > 0. This process is repeated for each position x s on the grid. For most locations, In 1Z is zero, but 
local maxima will occur in the vicinity of one or more events. Note that while the method is in some sense 
binned due to the discrete array of grid points, the spacing is chosen small enough so that errors introduced by 
binning are insignificant. 



Figure 1. Likelihood ratio In TZ, maximized with respect to n s , as a function of right ascension and declination for the 
combined set of HiRes events above 1.0 x 10 19 eV and AGASA events above 4.0 x 10 19 eV. 

For the signal probability density function (Q) of the HiRes events, we use a two-dimensional Gaussian func- 
tion whose width is chosen such that 68 % of the probability density function falls within an opening angle 
0.6°. Note that for a two-dimensional Gaussian distribution the opening angle 9 = 1.515 a encloses 68 % of 
the distribution. Since the dependence on energy is weak, we use the same value, o — 0.4°, for every HiRes 
stereo event. The procedure to obtain the background expectation is described in Q. 

3. Results 

We apply this method to the 271 HiRes stereo events with energies above 1.0 x 10 19 eV. The point with the 
largest In 1Z is at right ascension a = 219.2° and declination 5 = 74.9°. The maximum likelihood ratio at this 
position is In 1Z = 8.45 for n s = 2.9. 

The statistical significance of the appearance of a "source" with a maximum likelihood ratio In 1Z can be 
evaluated using simulated random data sets. The full likelihood analysis is performed for a large number of 
random data sets with the same number of events and the same underlying exposure as the original data set, 
but isotropic arrival directions. The chance probability for the "source" to appear is then given by the fraction 
of random data sets which have at least one location causing the maximum In 7Z to be equal or larger than 
8.45, the value of the maximum in the real data. For the HiRes stereo data above 1.0 x 10 19 eV, the chance 
probability of the source hypothesis is of the order of 28 %. Consequently, there is no statistically significant 
evidence for clustering consistent with a point source in the data set. 

We now include the published AGASA data set above 4.0 x 10 19 eV in this analysis. For AGASA, we ap- 
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proximate the probability density by the sum of two Gaussian functions chosen such that the 68 % and 90 % 
opening angle given in |7| is correctly reproduced. Details can be found in Q][S). 

Fig-E shows the result of the analysis. At each a and 6, the likelihood ratio is shown for the number of 
source events n s which maximizes In 1Z. The point with the largest In 1Z is at right ascension a = 169.0° 
and declination S = 56.2°, which is indeed the position of the point source claimed in |3|. The maximum 
likelihood ratio at this position is In 1Z = 8.38 for n s = 3.6. As before, we evaluate the chance probability for 
the appearance of a source with maximum In 1Z = 8.38 or higher in this data set by analyzing a large number 
of simulated isotropic data set, now containing 57 AGASA events in addition to the HiRes stereo events. The 
chance probability is 43 %, so again, no evidence for clustering consistent with a point source is found. 

The HiRes experiment is scheduled to take data until March 2006, and the search for cosmic ray point sources 
at the highest energies will continue when more data becomes available. 
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